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Introduction 

The problem of determining bearing loads and 
stability properties of rotating machines such as the 
turbopumps used in high performance rocket engines 
1 ike the Space Shuttle Main Engine (SSME) is compl ex . 
Very high speeds are attained with significant fl uid 
flows. As a consequence, bearing loads are poten- 
tially high with subsynchronous whirl ing likely. 
Typically, model s used to analyze such systems are 
very compl icated and nearly Impossible to use for 
gaining insight into the basic phenomena involved. 
Linear model s containing large numbers of degrees of 
freedom have been devel oped and applied to the analy- 
sis with mixed success. A significant nonl i neari ty 
is Ignored by these model s . The bearings typically 
have cl earances of the order of .0005 "-.0025". Since 
these machines are bal anced to very high precision, 
the eccentricity of the rotor, i.e. the distance be- 
tween the rotor center of mass and its geometric axis 
is of the same order or smaller in magnitude. Thus, 
bearing cl earances or deadbands as they are more typ- 
ically called, significantly affect the dynamics of 
these systems and must be taken into account. Taking 
this nonl i neari ty into account makes the analysis of 
the dynamics much more difficul t. It is very de- 
sirable to have a simpl ified model of a turbopump 
which retains the significant driving forces known to 
be present but readily lends itself to analysis. 
Such a model is available and is usually referred to 
as the Jeffcott model . We have modified this model 
by adding deadband effects along with fluid seal 
forces as currently understood. Further, we have 
rewritten the equations of motion for the model in 
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polar coordinates. This formulation is more naturally 
suited to the symmetry of the probl em because the 

whirl orbits tend to be circul ar. 

In addition to seal forces and deadbands, we 
have added a constant side force to the model to 
account for the 1 ikely mi sal ignments between bearings 
and seal s and al so to account for hydrodynamic forces 
resul ting from pumping fluids which may not be per- 
fectly bal anced due to si ight imperfections in the 

internal geometry of the pump. The side force and 
deadband effects, working together, significantly 
affect the stability properties of the system in an 
interesting way. Stabil ity may be enhanced under 
proper combinations but is only 1 ocal stabil ity in 

that it is possible to drive the system into insta- 
bil ity by impulsive disturbances or large rotor 

imbal ances. 

The Jeffcott rotor is cl oser to real i ty than it 
may appear to the casual observer. Periodic synchro- 
nous or nonsynchronous orbiting motions of the rotor, 
referred to as whirls, are normal ly the motions of 
the system exhibited. Such an orbital motion can be 
described by a planar model . Thus, values for the 
effective mass, stiffness, deadband and seal coef- 
ficients can be found which will approximate the 
behavior of the more compl ex models. While exact 
frequencies of critical speeds and stabil ity boun- 
daries cannot be inferred from Jeffcott model s, very 
good qual itative behavior can be investigated with 
these model s and ref i ned by higher f i del i ty hybrid 
simul ations. For this reason, we consider the 
augmented Jeffcott model as the model of choice for 
devel oping an understanding of rotor whirl 
phenomena. 

Force Models & System Equations of Motion 


The assumed geometry for the derivation of the 
equations of motion of the simple rotor model is 
depicted in Figure 1. The vector r_ is the displ ace- 
ment of the rotor center from its equil ibrium posi- 
tion (rotor at rest) . The angle <(> is the angle made 
by r with the horizontal axis and is referred to as 
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Figure 2. Vector force diagram. 



the whirl angle. The rotor eccentricity is repre- 
sented by the vector _e, the magnitude of which is 
constant. The shaft speed is denoted by oo and is 

assumed constant for the analysis. The equations of 
motion are derived in polar coordinates. The unit 
vectors _u rs ua, and jj* in the figure indicate the 
polar coordinate reference frame. 

The forces which must be considered in the for- 
mulation of the equations of motion include bearing 
forces, seal forces, imbalance forces, and side 
forces. Figure 2 is a vector force diagram which in- 
dicates the direction in which each of the various 
forces acts. The turbopump rotors are maintained in 
position by bearing forces. These bearing forces are 
generated by a rather compl ex interaction involving 
bending forces of the rotor shaft, the deformation of 
the bearing ball s or roll ers, the motion and defor- 
mation of the bearing races, the bearing retainers, 
the bearing carriers, etc. For our purposes, we as- 
sume that the bearing acts as a 1 inear spring. How- 
ever, clearances between bearing races and carriers 
or shafts allow some small region of free motion of 
the rotor shaft relative to its housing. For sim- 
pl icity, we ideal ize the bearing ball s or rol 1 ers as 
a uniform annul ar ring separating the rotor shaft and 
housing. The bearing force curve is ideal ized for 

the analysis as 

Kb (r - 9 u r ) lx I > 9 

£b = - — (i) 

1.0 |r|<g 

where Kr is the bearing stiffness and g is the dead- 
band. If the magnitude of £ is less than g, then the 
bearing forces are assumed to be zero. 

The fluid bei ng pumped reacts upon the rotor 
with forces that are dependent upon rotor position 
and vel ocity and can be represented by 1 inear model s 
for small displacements. 1 The seal s prevent the high 
pressure fl uid from leaking away and al so generate 
forces on the rotor which can be model ed 1 i nearly. 
The assumed form representing these forces is given 
by 
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fseal = -Csi - K s£ + Qs x I + C Q £x x I 


( 2 ) 


where K s is the seal stiffness, C s is the seal damp- 
ing coefficient, Q s is the cross coupling stiffness, 
and Cn is the cross coupling damping. These forces 
have the potential to drive whirl instability. 

The force due to the mass eccentricity is a 
rotating force whose magnitude varies as the square 
of the rotor speed and is directed toward the rotor 
center of mass. This force is potentially destruc- 
tive and must be minimized by stringent balancing of 
the turbopump rotors. The form of this force is 

Fg = -m 0 )^ _e (3) 

with m representing the mass of the rotor. The 
values of the parameters in the equations above are 
chosen such that the system model is representative 
of the SSME high-pressure oxygen turbopump. These 
values are: Kg = 10 6 lbs/in; K$ « 2.0 x 10 5 lbs/in; 

C$ = 200 lbs- sec/in; Cn = 40 lbs- sec/in; Q$ = C s cu/2 
lbs/in; and m = 0.20422 lbs-sec 2 /in. Five deadband 
values are considered and are 0.5, 1.0, 1.5, 2.0, and 
2.5 mils. Two values of rotor eccentricity are also 
used, they are 0.1 and 0.2 mils. Side force values 
range from 600 to 1200 pounds in the investigation. 

The force equation for the system may be wri tten as 

m(r + e) = I Fi , (4) 

_ _ 1 - 


with the summation representing the forces due to the 
seals, bearings, and the side force. The vector 
equation for the system is, therefore 

mr = -Kg (_r -guj) -K s _r - C s jr + 

n 

Qs jfx x X +C Q iJx x i + £s + m “’ £• (5) 

Performing the vector cross products and derivatives 
Indicated in equation 5 above, the fol 1 owi ng second 
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order differential equations which describe the 
system are obtained. 


" k B . . Ks C s # Cq # F s 

r = (r-g) r r — - r<ji + — cos<j> + 

m m m m m 

o 

rf + a) e cos (wt-<|>) 

.. _ Qs Cs . + Cq ^ F _± sin ^ 

m m m r m r 

2 

cj e 

+ sin (out - <j>) 

r 


(6) 


(7) 


Limit Cycle Analysis 

Analysis of the system is more readily carried 
out when the equations of motion are cast in state 
variable form. 


Let 


PI = r p 2 = r P3 


Then 
PI = P2 


p4 = <f> (8 ) 

(9) 


-(K B + K s ) C s Cq 

P2 = Pi + KB9 P2 Pi P4 

m mm 


r S ? 9 

+ COS P3 + oj e cos (out - P3) + pi + P4 
m 

P3 * P4 


( 10 ) 


( 11 ) 
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Cq P2 F s 


sin p 3 


Qs C s 

P4 P4 + 

mm m pi m pi 

( 12 ) 


2 

on e 

+ sin (out- P3 ) 

PI 


2 P2P4 
Pi 


In this form, the system is amenable to solution by 
numerical methods. Three different types of orbits 
are observed when the system of equation is solved 
for different values of deadband, side force, and 
rotor eccentricity. These are referred to as A- type, 
B-type, and C-type motion. 2 Illustrations of the 
three are shown in Figures 3 through 5. A- type 
motion is periodic and does not encircle the origin. 
B-type motion is nonperiodic and somewhat random in 
nature. C-type is periodic and does surround the 

origin. The A- type and C-type motions are 1 imit 

cycles. 

To characterize the limit cycle motions present 
in the rotor system, an algorithm has been developed 
which will converge to a set of initial conditions 
for the four system states which, when input into a 
simul ation of this system, will cause the system to 
immediately exhibit the limit cycle behavior. The 
al gori thm is based on the fact that the function for 
which the limit cycle initial conditions are sought 
is periodic. That is, the orbit comes back around to 
the same poi nt once each cycle. The idea is to 
determine that such a point exists and the values of 
the system states which sati sfy this condition. 

Given the state equations which describe the 
system, a sol ution to the states may be obtained 
through integration. The mathematical statement of 
the probl em is: 

p_ = (£,t) (13) 

the solution to which is 
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Figure 5. C-tyoe rotor notion. 
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(14) 


t 

_p(t) = £ 0 + / _£. (£,x)dx 

0 

where £ 0 is some initial state. It is desired to 
determine the j> 0 such that 

j>(t) = _p 0 for t = T (15) 

2ir 

where T = — , (16) 

0 ) 

the peri od of the function. In other words, we wish 
to determine j? 0 so that the integral in Equation 14 
is zero. The problem may be restated as 

_p(T ) = J> 0 + £(£) • (17) 

If £(£) can be driven to zero, then £(T) = £ 0 . 

The function £(£) may be approximated to 1st order by 

3g 

g(£) = g (£ 0 ) + — 

9 £ 

where A£ is some incremental change in the state vec- 
tor p. This is the quantity to be determined. It 

will've added to the original state vector. Because 

we wish £(£) to be zero. Equation 17 is rewritten as 

£ ■ £ (£ 0 ) + J • (£ - £ 0 ) (19) 

where 

A£ = £ - £ 0 (20) 

and J is the Jacobian of £(£). The sol ution for Ap 
is 
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A_P = J -1 (-l(p 0 )) . 


( 21 ) 


A new set of initial states is formed as 

Po new = Pool d + A P * (22) 

For orbits which are C-type, another modification to 
the algorithm is required. Because the whirl angle, 
<l> , is not periodic but increasing with time, this 
method left unmodified will not converge to a solu- 
tion. To force <j> to appear to be periodic, the value 
2tt is subtracted from the <j> component of the state 
vector value at time t = T. Thi s procedure will, in 
fact, allow the al gori thm to converge to a sol ution 
to the C-type orbit initial conditions. 

The plot of the 1 imit cycle shown in Figure 3 is 
obtained usi ng the initial conditions obtained with 
the algorithm described above. Figure 6 is a pi ot of 
that same 1 imit cycle including the transients that 
occur when a simulation of the system is executed 
with the initial conditions set at zero. Likewise, 
the plot in Figure 5 is obtained using 1 imit cycle 
initial conditions obtained with the algorithm. 
Figure 7 is a plot of the same orbit obtained from 
zero initial conditions. 

Stabil ity Analysis 

The approach taken in the determination of system 
stabil ity is to examine the equations of motion, 
omitting the imbal ance force terms. This procedure 
greatly simpl ifies the analysis. Our studies indi- 
cate that side forces and deadbands are infl uential 
on local stabil ity. Under the infl uence of a side 
force, the rotor shi fts to a position of equil ibrium 
when no imbal ance is present. The addition of an 
imbal ance, in general , causes the rotor to whirl 
about the equil ibrium point. The position of the 
equil ibrium point is dependent upon the magnitude of 
the side force, the deadband, and the stiffness coef- 
ficients. The position of the equil ibrium, with 
respect to the deadband, determines the type of orbit 
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in which the rotor will whirl. Generally, if the 
equilibrium point is outside the deadband, the rotor 
exhibits A- type motion with C-type motion resulting 
if the equil ibrium point is inside the deadband. 

The two parameters defining the equil ibrium point 
are rg, the magnitude of the displ acement of the 
rotor center from the rest position, and <|>g, the 
angle made by rg with the horizontal axis. If the 
rotor is in equilibrium, then the forces in the 
radial and transverse directions must both sum to 
zero. Because they have been bal anced by the side 
force, force terms arising due to the spin of the 
shaft have no infl uence on the equil ibrium point. 
The two equations defining the equil ibrium point are, 
therefore, 

Kb (r 0 - g) + K s r 0 = F s cos <j> 0 (23) 

Qs r o " F s sin <j> 0 (24) 


where F s is used to denote side force magnitude. If 
we define F r to be the side force in the radial 
direction and F^ to be the side force in the 
tangential direction, then the fol 1 owi ng expression 
is true. 


F s 2 = F r 2 + F/ (25) 

Using equations 23 and 24 with equation 25 makes 
possible the sol ution for rg which is 


no 


(K B +K s )K B g ± ^F s 2 [(K B+Ks ) 2 +Q s 2 ] - Q S 2 K B V 


(K B + K s ) 2 + Q s 2 


(26) 


Note that when there is no deadband, any value of 
side force will resul t in a positive sol ution for rg. 
However, in the presence of a deadband, the magnitude 
of the side force must be sufficient to push the 
rotor out to the bearing or el se the bearing stiff- 
ness plays no role in the determination of rg. The 
minimum side force required to overcome the seal 


356 



stiffness is denoted F$min ancl 1S defined by the 
expression 

F$MIN " V (K s 2 + Q s 2 )g 2 ' • (27) 

For any value of side force less than F$min> the 
value of rg is the positive solution to Equation 26 
with set to zero. For cases in which the side 
force is greater than F$min and two positive solu- 
tions result, then rg is equal to the larger of the 
two val ues because it will 1 ie outside the deadband 
due to the side force being sufficient to move the 
rotor to such a position. The angle <f>o is easily 
obtained once the sol ution for rg is determined. 

With the equil ibrium points well in hand, we may 
now proceed. Recall the vector equation which 
describes the system. It is repeated here, with the 
imbal ance term omitted. 

mr = -Kb( r - g) u (r - g) e r -K s re r + Qge* x re r 


- C s jr + Cq< 2 * x _r (28) 

We make the fol 1 owi ng definitions for _r and _e r , 

T_ = £ 0 + _S (29) 

r = r 0 + 6y (30) 

e r = e^ 6 e r ( 31 ) 


where r 0 is the equil ibrium position vector, ej Q is 
the unit vector in the direction of _r 0 , and 6_ and 6_er 
are the perturbations associated with £ and £ r> 

respectively. The radial component of jr with its 
perturbation is Equation 30. Another way to express 
£ r is 

lo + i lo A lo lo 

_ep = = + - • 6_ 

I lo + £l I Ho I I lo i I r 0 | 2 

+ ... (32) 
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which is equivalent to: 


$2 

Oy* = ■“*- Q( j> « ( 33 ) 

r 0 “ 

We now examine the nonlinear deadband term in 
Equation 27. In the small and to a first order 
approximation, 

-Kb (r - g) u (r - g) e r 

(r 0 - g) (34) 

= -Kerqero + K B9£r 0 ■ K B 6 y e r 0 ~ K B 6 2 £<j> 

r o 

with Sy and 6 Z being the small perturbations about 
the equil ibrium. It fol 1 ows , therefore, that in 
terms of the perturbation variables, the system may 
be expressed in the fol 1 owi ng form: 



9 

Where g 0 = — . 
ro 


The effects of side forces are now inherent in the 
formulation. Stabil ity may be assessed through exam- 
ination of equation 35. State assignments for the 
perturbation variables are given bel ow. 
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X 1 = 6 y x 2 6 y X 3 = 6 Z X4 = 6 2 (36) 

Rewriting 34 in state variable format yields the 
following differential equations. 


X 1 = *2 

. K C s Q$ Cq 

X 2 = X X x 2 - x 3 - X 4 

m m m m 

« 

x 3 = X 4 

. Q$ Cq k s + Kb (1 - go ) 

*4 = Xl + x 2 x 3 

mm m 


(37) 

(38) 

(39) 

(40) 


Cs 

X4 

m 


The sum Kg + k s has been repl aced by K for sim- 
plicity's sake. The system matrix is formed as 
before from which the characteristic equation is 
derived by solving the determinant of [si - A]. 


P(s) = 



Cs 

+ — s + K + 

m 




Cs 

+ — s + K - 
m 



( 41 ) 


The stabil ity boundaries establ i shed are pi otted in 
Figure 8. Several interesting facts are observed 
when this figure is examined. The stabil ity boundary 
for a deadband of zero is a constant 4848 radians/ 
second, the dashed curve in the figure. This is the 
same stability boundary which may be establ i shed for 
the simpl est form of the system. 3 This frequency is 
considered to be the global stabil ity boundary. That 
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is, the system is globally unstable when run at fre- 
quencies higher than this value. Local stabil ity, or 
stabil ity in the small, may be enhanced by other fac- 
tors. 

Notice that all five of the non-zero deadband stabi- 
1 ity curves are very similar in their general charac- 
ter. The maxima appear at approximately the same 
spin frequency as do their minima. We have shown 
that these curves do, indeed, collapse into a si ngl e 
curve when the system is nondimensional ized. To ac- 
complish this, the units of displ acement, force, and 
time are modified in such a way that the system 
parameters become unitl ess. Di spl acements are ex- 
pressed in units of g, the deadband. Forces may be 
expressed as the product of the seal sti ffness and 
deadband. Time is exp ressed in term s of the system 
natural frequency is Vm/(K S + Kg) »’ Figure 9 is a 
plot of the stabil ity boundary for the normalized 
system. With this curve and the given conversion 
factors, one may determine the stabil ity boundary for 
any deadband value. 

Bearing Loads Considerations 

One of the major points of the study is to deter- 
mine the effects of the system parameters on bearing 
loads. If these loads become too large, the effects 
are detrimental . Initially, we will look into this 
matter using a system with no side forces present, 
but with a rotor imbal ance. We will then take into 
consideration the added effects of a side force pre- 
sent in the system. 

PI otted in Figure 10 are the bearing loads which 
resul t when the rotor eccentricity is 0.1 mil s. The 
deadband range is from 0.0 to 0.2 mil s. The general 
behavior is that the smal 1 er deadband produces the 
1 argest bearing load. This makes sense because the 
seal forces must be overcome before there is any 
i nteraction between the rotor and the bearings. The 
more distance between the rotor and the bearing there 
is, the more effect the seal forces have. It has 
been shown that if the rotor eccentricity is doubled, 
the bearing loads al so double if all other parameters 
are left unchanged. 3 
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Figure 10. Bearing force curves. 



Figure 11. Bearing force curves. 
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Bearing load analysis is performed for two non- 
zero side force values. A value of side force is 
chosen so that it is al ways greater than F$mjn for 
any of the five deadbands considered up to a fre- 
quency of 5000 radians per second. Figure 11 is a 
plot of the bearing loads for the deadbands of 0.5 to 
2.5 mils for the side force of 1350 pounds. The max- 
imum bearing loads occur at the system natural fre- 
quency of approximately 2424 radians/second, the 
smal 1 est deadband producing the largest load. The 
load curves are pi otted only up to a shaft spin fre- 
quency of 4000 radians because the system becomes 
unstable for frequencies higher than that. The pres- 
ence of the rotor eccentricity of 0.2 mil s is respon- 
sible for the unstable behavior. 

A simil ar family of curves is produced when the 
side force is increased to twice that of the "minimum" 
side force, or, 2700 pounds. The plots of these 
bearing loads are given in Figure 12. The same gen- 
eral behavior is exhibited as before. Instabil ity 
occurs somewhat sooner, at 3400 radians/ second. The 
1 oads are much greater, as well . 

The curves presented were generated using the 
simulation to determine the maximum rotor displ ace- 
ment, after steady-state is achieved. Having this 
val ue, it is a simpl e matter to compute the bearing 
load. Equation 43 is used. 

BFMAX = MrMAX - 9 ) (43) 

Conclusions 


In the previous sections we have discussed the 
model ing and analysis of the problem. It is now time 
to review and summarize our resul ts. 

1. Observed 3 motion types cal 1 ed A, B, C; 

A - Periodic but does not enclose origin, may 
i ncl ude higher harmonics; 

B - Nonperiodic; 

C - Periodic enclosing origin, synchronous or 
nonsynchronous; 
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2. Limit Cycle Algorithm developed and employed, 
both A & C types observed. 

3. Deadband does not affect stab i 1 i ty- i n- the- 

1 arge. 

4. Stability-in- the- small are affected (enhanced) 
by deadband and si deforce. 

5. Bearing loads are 1 argest for C-type motion. 

6. Side force acting in concert with deadband 
effects may either increase or decrease bearing 
loads. 

7. Bearing loads in a stable pump are determined 
primarily by rotor imbal ance and side forces. 

These resul ts are quite significant in our under- 
standing of the effects of bearing deadbands. 
Harmonics of snychronous and nonsynchronous oscilla- 
tions have been observed. This is cl early a nonl in- 
ear effect. Stable limit cycle whirl s have been 
observed occurring at synchronous and nonsynchronous 
rotor speeds in our resul ts. 

The limit cycle algorithm that we have developed 
can be general ized to more complex turbopump model s 
with more degrees of freedom. It will be useful for 
loads analysis with nonl i near forces for rotor dyna- 
mics and other appl ications. It is capable of con- 
verging to periodic motions (sol utions) which 
generally resul t in the highest load- producing con- 
ditions. 

Since stabil ity-in-the-large is ul timately deter- 
mined by behavior at extremely large ampl itudes of 
motion, deadband effects become negligible. Thus, 
1 i near model s remain adequate for analysis of global 
stabil ity properties. Stabi 1 1 ty- i n- the- smal 1 is 
significantly al tered by the nonl inear effects of 
deadbands. We have shown that sideforces can signi- 
ficantly enhance stabil ity provided imbal ance offsets 
and/or impul sive disturbances do not cause signifi- 
cant displ acement from the equilibrium position of 
the rotor. 

Bearing loads have been shown to be significantly 
modified by deadband effects. Critical speeds are 
altered. Loads may increase or decrease. The shape 
of the critical response curve is al tered with higher 
1 oading at lower frequencies due to the deadband. 
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These results have been obtained using a relati- 
vely simple 2 degree- of- freedom model. This may lead 
one to believe the results are not applicable to real 
machines. This Is not the case, however, and Indeed 
one can argue and demonstrate with more sophisticated 
model s that these effects are real . Si nee rotor 
responses are most often periodic, such motions can 
be described adequately by an effective mass 
responding to effective stiffnesses and deadbands, 

l.e., a 2-dimenslonal model. Thus, our results are 
at least qualitatively valid for the description of 
turbo-pump motions. 
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